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We propose a simple realization of the three-dimensional (3D) Weyl semimetal phase, utilizing 
a multilayer structure, composed of identical thin films of a magnetically-doped 3D topological 
insulator (TI), separated by ordinary- insulator spacer layers. We show that the phase diagram 
of this system contains a Weyl semimetal phase of the simplest possible kind, with only two Dirac 
nodes of opposite chirality, separated in momentum space, in its bandstructure. This Weyl semimetal 
has a finite anomalous Hall conductivity, chiral edge states, and occurs as an intermediate phase 
between an ordinary insulator and a 3D quantum anomalous Hall insulator. We find that the Weyl 
semimetal has a nonzero DC conductivity at zero temperature, but Drude weight vanishing as T 2 , 
and is thus an unusual metallic phase, characterized by a finite anomalous Hall conductivity and 
topologically-protected edge states. 



The recent discovery of time-reversal (TR) invariant 
topological insulators (TI) pQ has led to a surge of in- 
terest in topological properties of the electronic band- 
structure of crystalline materials. TIs exhibit a bulk gap, 
but gapless surface states, whose gaplessness is protected 
by topology. Remarkably, recent work has demonstrated 
that such a surface-bulk correspondence can also obtain 
even when the bulk is gapless, by virtue of point touch- 
ings of non-degenerate conduction and valence bands [5] . 
Such accidental point touchings have been known to exist 
since the earliest days of the theory of solids [3] , but their 
topological properties have only been appreciated much 
more recently [2JI3HE], and concrete materials, where they 
may be found, have been proposed [SMS]. Non-trivial and 
robust band touching requires either broken TR or inver- 
sion (I) 0E1E], in which case the touching points acquire 
topological character and thus give rise to stable phases 
of matter. The bandstructure near these points can be 
described by a massless two-component Dirac or Weyl 
Hamiltonian: 



H = zkvpcr ■ k, 



(1) 



where k is the crystal momentum in the first Brillouin 
zone (BZ), expanded near the band-touching point, cr is 
the triplet of Pauli matrices and the sign in front corre- 
sponds to two different possible chiralities, characterizing 
the point. Such Weyl fermions have been studied exten- 
sively in high-energy physics, in particular as a descrip- 
tion of neutrinos [J], and may be viewed as topological 
defects (hedgehogs) in momentum space [1] . Any pertur- 
bation of Eq.Q only shifts the degeneracy point in en- 
ergy or momentum, but does not remove it: an isolated 
Weyl fcrmion in this sense possesses an absolute topo- 
logical stability (this is in contrast to 2D massless Dirac 
fermions in graphene, where inversion symmetry of the 
honeycomb lattice is essential for their stability). Very 
general considerations show that Dirac degeneracy points 
can only occur in pairs of opposite chirality [7] and can 
thus be eliminated by pairwise annihilation. When the 



TR or I symmetry is broken, however, the Weyl fermions 
are separated in momentum space and thus, assuming 
translational symmetry remains intact, are still topolog- 
ically stable. 

Ref. [5] has proposed a possible realization of a Weyl 
semimetal with 24 Dirac nodes in iridium pyrochlores, 
which are strongly-correlated magnetic materials (a dif- 
ferent scenario for this material was proposed in Ref. [B]). 
The purpose of this work is to propose a much sim- 
pler realization of the Weyl semimetal, not relying on 
strong correlations in a rather complex material. The 
Weyl semimetal we propose also possesses only two Dirac 
nodes, the smallest possible number, and thus is in a 
sense the most elemental realization of this phase of mat- 
ter. 

The material we propose is a multilayer heterostruc- 
ture, consisting of alternating layers of a 3D TI material, 
such as Bi2Se3, and an ordinary insulator, which serves 
as spacer material between the neighboring TI layers, 
as shown in Fig. [T] The ability to grow ultrathin high- 
quality films of Bi2Se3 has been clearly demonstrated in 
recent experiments [9] . It is thus quite realistic to expect 
that a multilayer heterostructure, consisting essentially 
of a stack of such thin films, can be fabricated using 
available technology. The Hamiltonian, describing this 
heterostructure, can be written as: 



H = [v f t z (z x cr) • k±Si,j + ma z Si_j + AsT x S lyJ 



(2) 



The first term in Eq.Q describes the two (top and bot- 
tom) surface states of an individual TI layer. We as- 
sume for simplicity that a TI material with a single two- 
dimensional (2D) Dirac node per surface BZ is employed. 
vf is the Fermi velocity, characterizing the surface Dirac 
fermion, which we take to be the same on the top and 
bottom surface of each layer, kj^ is the momentum in 
the 2D surface BZ (we use Ti = 1 units), cr is the triplet 



FIG. 1. Schematic drawing of the proposed multilayer struc- 
ture. Unhashed layers are the TI layers, while hashed layers 
are the ordinary insulator spacers. Arrow in each TI layer 
shows the magnetization direction. Only three periods of the 
superlattice are shown in the figure, 20-30 unit cells can per- 
haps be grown realistically. 

of Pauli matrices acting on the real spin degree of free- 
dom and r are Pauli matrices acting on the which sur- 
face pseudospin degree of freedom. The indices i,j label 
distinct TI layers. The second term describes exchange 
spin splitting of the surface states, which can be induced, 
for example, by doping each TI layer with magnetic im- 
purities, as has been recently demonstrated experimen- 
tally [TU]. The remaining terms in Eq.Q describe tun- 
neling between top and bottom surfaces within the same 
TI layer (the term, proportional to Ag), and between 
top and bottom surfaces of neighboring TI layers (terms, 
proportional to Ad). Longer-range tunneling is assumed 
to be negligible. We will regard to and As,d as tunable 
parameters and study the phase diagram of Eq.Q as a 
function of these parameters. 

Let us initially assume that the spin splitting is ab- 
sent, i.e. set to = 0. Diagonalizing Eq.Q one finds the 
following band dispersion: 

el ± = v 2 F (k 2 x + k 2 y ) + A 2 (k z ), (3) 

where A(k z ) = y/A 2 s + A 2 D + 2A S A D cos(k z d), and d is 
the superlattice period (i.e. TI layer plus spacer layer 
thickness) in the growth (z) direction. This bandstruc- 
ture is fully gapped when |Ag| ^ |Arj|, but contains 
Dirac nodes when As/Ap = ±1. The nodes are lo- 
cated at k z — ir/d when As/ An = 1 and at k z = 
when As/ An = — 1 (k x = k y = always). While both 
cases are possible, we will assume the former for concrete- 
ness and will take both tunneling matrix elements to be 
positive (this choice does not affect any of our results). 
Expanding the band dispersion near the Dirac point at 
k x = k y — 0, k z — ir/d to leading order in the momentum 
one obtains: 

el ± =v 2 Akl + k 2 y )+v 2 F k 2 z , (4) 

where v F — d^/AsAp. The momentum-space Hamilto- 
nian near the Dirac node has the form: 

n(k)^v F T z {zxa)-k + d F T y k z , (5) 

which can be brought to a block-diagonal form, explic- 
itly revealing a pair of two-component Weyl fermions 



(a) m=0 (b) m*0 

FIG. 2. (Color online) Phase diagrams for (a) m = and 
(b) m 7^ 0. In (a), the red line represents the phase bound- 
ary between topological insulator (TI) and ordinary insulator 
(Ins). In (b), due to TR symmetry breaking, the distinction 
between topological and ordinary insulators is moot, so the TI 
in (a) has been converted to Ins. QAH denotes the quantum 
anomalous Hall phase. 

with opposite chirality, by a n/2 rotation around the x- 
axis in the pseudospin space. Alternatively, in total this 
is a conventional 4-component massless Dirac fermion. 
As discussed above, since the two Weyl fermions are lo- 
cated at the same point in momentum space, they are 
topologically unstable. Any perturbation, for example 
any deviation of the ratio As/Ap from unity, immedi- 
ately eliminates the degenerate Dirac node and produces 
a fully gapped spectrum. With m = 0, the massless 
Dirac point can be understood [S] as a critical point be- 
tween topological (An > As) and ordinary (An < As) 
insulators with both inversion and time-reversal symme- 
try preserved (see Fig. 2a). To produce a topologically 
stable phase with 3D Dirac nodes, the nodes have to be 
separated in momentum space. As mentioned above, this 
can generally be accomplished by breaking either TR or 
I symmetries and there are in principle many ways to do 
this. Here we will focus on one particular way, which is 
perhaps the simplest from the point of view of a practi- 
cal realization. Namely, as already mentioned above, we 
will assume that each TI layer is doped with magnetic 
impurities, producing a ferromagnetically-ordered state 
within each layer, with magnetization along the growth 
direction of the heterostructure. This leads to spin split- 
ting of the surface states of magnitude to, described by 
the second term in Eq.([2]). The band dispersion is now 
given by: 

4± =v F (kl + k 2 y ) + [m±A(k z )] 2 . (6) 

This dispersion has two nondegenerate Dirac nodes, sepa- 
rated along the z-axis in momentum space, with locations 
given by k z = ir/d ± fen, where: 

k = - arccos [l - (to 2 - (A s - A D f)/2A s A D ] . (7) 

The nodes exist as long as: 

TO 2 ! = (A s - An) 2 < TO 2 < TO 2 2 = (A s + An) 2 . (8) 
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Thus, as discussed above, splitting of the degenerate 
Dirac node in momentum space produces a stable Weyl 
semimetal phase, existing in a finite region of the phase 
diagram (Fig. 2b). The Weyl semimetal occurs as an in- 
termediate phase between an ordinary insulator (to 2 < 
to 2 x ) and a 3D quantum anomalous Hall (QAH) insula- 
tor [TT] with quantized Hall conductivity, equal to e 2 /h 
per TI layer (to 2 > to 2 2 ). As to is increased from zero, a 
degenerate Dirac point at k x = k y = 0, k z — -k / d appears 
at the lower critical value (rn c i) of the spin splitting. 
The Dirac node is split along the z-axis into two non- 
degenerate nodes in the Weyl semimetal phase, with the 
splitting increasing monotonically with the magnitude of 
to. At the upper critical value (m C 2), the two nondegen- 
erate nodes meet and annihilate at the center of the BZ 
k x = k y = k z = 0, giving rise to the fully gapped QAH 
insulator. 

We will now show that the Weyl semimetal phase we 
have found is characterized by a nonzero, but nonquan- 
tized anomalous Hall conductivity, proportional to the 
magnitude of the separation of the Dirac nodes in mo- 
mentum space, thus smoothly evolving from zero in the 
ordinary insulator phase to e 2 /h per layer in the 3D QAH 
insulator phase. 

The fact that a 3D Weyl semimetal generally has a 
finite Hall conductivity is known, and was pointed out 
e.g. in Ref. [T2]. For the particular realization of a Weyl 
semimetal, proposed in [2], the Hall conductivity van- 
ishes due to the cubic symmetry of the crystal structure 
of the pyrochlore iridate materials. If the cubic symme- 
try is broken, e.g. by applying a uniaxial pressure, as 
proposed in [BJ, the Hall conductivity becomes nonzero, 
and is proportional, in the case discussed in |BJ, to the 
applied pressure. In our case, there is a natural pre- 
ferred axis, i.e. the growth direction of the heterostruc- 
ture (which coincides with the magnetization direction) 
and the Hall conductivity in the Weyl semimetal phase 
is automatically nonzero. The simplest and most physi- 
cally transparent way to obtain this result is to view the 
3D bandstructure as a set of independent 2D bandstruc- 
tures at fixed k z . We begin with the momentum-space 
Hamiltonian: 

HQs) = v f t z (z x cr) • k + ma z + A(k z ), (9) 

where A = A s t x + \ (A D T+e ik * d + h.c). This is sim- 
plified by the canonical transformation: 

a ± ^r z a ± , T ± ^a z T ± . (10) 

After this transformation, the Hamiltonian becomes: 

HQs) = v F k y o x - v F k x a y + [to + A(k z )]a z . (11) 

Now A is a constant of motion, and may be replaced 
by its eigenvalues A(k z ) = ±A(k z ). For each of these 




FIG. 3. (Color online) A plot of a xy in units of e 2 /hd for 
a multilayer consisting of 30 TI layers. Ad is taken to be 
equal to 0.8As. The Hall conductivity is zero in the ordinary 
insulator phase (m < Ag - Ad) and reaches the maximum 
value of e 2 /hd in the QAH phase (m > As + Ad). Plateaus 
correspond to the wavevector ko being in an interval between 
the neighboring finite-size-quantized k z values. 



vanishes when k z = n/d ± fc , corresponding to the two 
Dirac nodes (M+ never vanishes). 

It is well-known that a change of sign of a 2D Dirac 
mass signals a quantum Hall transition, at which the 
quantized 2D Hall conductivity a 2 ® jumps by e 2 /h |13j . 
The absolute value of the Hall conductance is not deter- 
mined by the above continuum model. Therefore the con- 
tribution to the total 3D Hall conductance of the states at 
fixed k z is equal to cr 2 y(k z ) = e 2 /h {n+Q(k — \k z —Tr/d\)}, 
where Q(x) is the Hcaviside step function, and n is an 
integer. Since when to vanishes TR symmetry demands 
that the Hall conductivity must also vanish, we can con- 
clude that n — 0, and hence: 



*l d dk 

un 'z 2D 

-TT/d 27r 



e 2 k 
irh 



(12) 



Thus the anomalous Hall conductivity in the Weyl 
semimetal is proportional to the separation of the Dirac 
nodes in momentum space. For a multilayer, consisting 
of a finite number of layers, as will be the case in ex- 
periment, a xy will exhibit plateaus as a function of to, 
when &o will fall in an interval between the neighboring 
quantized k z values, as shown in Fig. [3j At the upper 
critical value of to — to c2 , 2k = 2ir/d, the two Dirac 
nodes annihilate each other at the center of the BZ, and 
the Hall conductivity reaches a quantized value per TI 
layer: 



e 



(13) 



cases, Eq. (Ill gives a 2D Dirac Hamiltonian for fixed 
k z , with a mass M± = to ± A(k z ). For to > m c i, Af_ 



which characterizes the to > m C 2 3D QAH insulator 
phase. 

The results of Ref. [5] imply the existence of "Fermi 
arcs" for the Weyl semimetal, in this case for any surface 
except the one normal to the z axis. In fact, this arc is 
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nothing but the set of edge states corresponding to the 
2D integer quantum Hall states for ir/d—ko < \k z \ < ir/d. 



These can be explicitly found from Eq. (Ill for e.g. the 



case of a surface at y = 0, modeled by a y-dependent spin- 
splitting field m(y). Outside the sample (y > 0) we take 
m {y ~~ * +oo) = 0, which realizes an ordinary insulator, 
while m > m c \ for y < 0, and consider the eigenstates for 
fixed k Xl k z . There are then special surface wavefunctions 
with A(k z ) — —A(k z ), which are eigenstates of (111: 



AnAk x ,k z ;y) = e J?VM»')-A(k,)]/»F| (r i, = (14) 

For n/d — ko < \k z \ < w/d and not otherwise, the ex- 
ponential above vanishes when y — > ±oo and the state 
is normalizable and localized to the surface with a local- 
ization length ~ vp /M-(k z ). The energy of this 
state is simply e sur f = vpk x , indeed identifying it as a 
chiral edge state. As m is increased to enter the QAH 
insulator phase, the arc extends across the full BZ and 
then the surface states can be alternatively viewed in a 
Wannier basis localized in individual TI layers, i.e. as 
conventional quantum Hall edge states for each layer. 

While the Weyl semimetal has topologically-protected 
edge states (topological protection follows from the sep- 
aration of the Dirac nodes in momentum space and from 
the fact the edge states are chiral), it is actually not a 
Hall insulator. At first sight, the low density of states 
g{e) ~ e 2 , seems to suggest a vanishing DC conduc- 
tivity at zero temperature [2j. However, a careful cal- 
culation reveals that, in the presence of disorder but 
not interactions, this is not the case. Using the stan- 
dard Kubo formula expression (or Boltzmann equation) 
with the Born-approximation impurity scattering rate 
l/r(e) = 2Trjg(e), where g(e) = e 2 /2ir 2 v F and 7 charac- 
terizes the strength of the impurity potential, the optical 
conductivity of a Weyl semimetal with isotropic Fermi 
velocity vf at temperature T, is given by: 



Re o-(ui) 



dx 



x 4 sech 2 (x) 



x 4 + (/i 3 4w/327r 2 7 T 2 ) 2 ' 



, (15) 

where we have restored explicit h for clarity. This gives a 
finite DC conductivity or>c ~ e vp/h^, which can be ex- 
pected to be large in a clean multilayer, but a Drude-like 
peak in the optical conductivity with weight, vanishing 
as T 2 . Thus with disorder (but neglecting interactions), 
the Weyl semimetal is not an insulator, but an unusual 
metal, characterized by a nonzero anomalous Hall con- 
ductivity and topologically-protected edge states. 

In conclusion, we have proposed a simple realization 
of a 3D Weyl semimetal phase in a multilayer structure, 
composed of a stack of thin layers of magnetically-doped 
3D TI material, separated by insulating spacers. We have 



shown that this material realizes the simplest possible 
type of Weyl semimetal, with only two Dirac nodes, sepa- 
rated along the growth direction of the heterostructure in 
momentum space. This Weyl semimetal is characterized 
by a nonzero anomalous Hall conductivity, proportional 
to the separation between the Dirac nodes, and by the 
existence of topologically stable chiral edge states. These 
edge states are, however, distinct from the ordinary quan- 
tum Hall edge states, since they exist not in the whole 
edge BZ, but in its finite subset, whose size is determined 
by the momentum-space separation of the Dirac nodes. 
Finally, we find that the Weyl semimetal has a finite DC 
conductivity at zero temperature, but Drude weight van- 
ishing as T 2 , and is thus an interesting metallic state, 
characterized by a nonzero anomalous Hall conductivity 
and topologically-protected edge states. Interesting open 
questions include the influence of Coulomb interactions 
on the properties of Weyl semimetals, in particular their 
transport properties. 
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